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ABSTRACT 


The problem of electrostatic wave generation by a return current driven 
by a small area electron beam during solar hard X-ray bursts is discussed. 

The marginal stability method (cf. Duijveman et al. 1981) is used to solve 
numerically the electron and ion heating equations for a prescribed beam 
current evolution. When ion— acoustic waves are considered, we find that the 
method appears satisfactory and, following an initial phase of Coulomb 
resistivity in which T g /T^ rises, predicts a rapid heating of substantial 
plasma volumes by anomalous ohmic dissipation. This hot plasma emits so much 
thermal bremsstrahlung that, contrary to previous expectations, the unstable 
beam-plasma system actually emits more hard X-rays than does the beam in the 
purely collisional thick target regime relevant to larger injection areas. 
Inclusion of ion-cyclotron waves results in ion-acoustic wave onset at lower 
Tg/Ti and a marginal stability treatment yields unphysical results. 
Specifically, negative resistivity occurs when the ion-acoustic instability 
is excited at an electron-to-ion temperature ratio of T /T- < 4»5, showing 
the marginal stability analysis of ion acoustic waves to be invalid in this 
regime. Discarding marginal stability and adopting a simple wave-energy 
equation, the time— dependent effect of ion— acoustic turbulence generated by 
the return current is investigated. 


1 Introduction 

It is commonly believed that electron beams, propagating downwards in the 
solar atmosphere, play a major role in the production, by collisional 
bremsstrahlung, of hard X-ray bursts during the impulsive phase of solar 
flares (see review by Brown and Smith 1980) and that they may also be 
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instrumental in flare atmospheric heating. In addition, it is recognised 
(e.g. Hoyng, Brown and van Beek 1976) that the large electron flux, demanded 
by the observed X-ray photon flux in such an interpretation, requires that a 
beam-neutralising return current be set up, such that: 

Hv V, 


n 


( 1 ) 


where 


^b 

n 


is the ratio of beam density to ambient plasma density, v, is the 


(sub-relativistic) beam velocity (a beam with single injection energy E Q , 
representing the mean energy of a real beam, is considered here) and v^ is 
the drift velocity of the ambient plasma electrons (which constitute the 
return current). If we consider, for simplicity, a model in which v d increases 
during the impulsive phase, due to the rising beam flux n b v b , threshold 
velocities for the generation of various microinstabilities may be reached and 
so unstable wave growth may take place in the atmosphere. The resultant 
anomalous resistivity will affect hard X-ray production in two ways: 


(a) The electric field required to drive the return current will increase 
and so reduce the lifetime of beam electrons. The nonthermal bremsstrahlung 
efficiency, which is already small, is therefore further reduced. 


(h) Enhanced plasma heating will take place and thereby affect the 
observed thermal radiation signature. 


The aim of this paper, then, is to describe preliminary results of 
calculations attempting to determine the effect of such return current 
instability on hard X-ray bursts during solar flares. 


2 Marginal Stability Analysis 

We consider a quasi-steady state consisting of a monoenergetic driving 
electron beam, with specified current density j b (t), current neutralising 
(jp = “ Jb) hot lifting electrons at temperature T 0> and hot stationary ions 
at temperature T^. While recognising the possibility of beam-return current 
interaction, we concentrate here on the return current electron-background ion 
instability. The quasi-linear relaxation of the electron beam is considered 
elsewhere in these proceedings by McClements et al. 

A typical large solar flare value is chosen for the peak total electron 
injection rate, ^(s" 1 ), but the beam area, A, is taken to be well below the 
upper limit set by hard X-ray images in order to ensure that unstable return 
current drift velocity thresholds are exceeded before this peak, viz. 

& Q - 10 36 s" 1 ; A = 10 16 cm 2 (2) 

(Brown and Melrose 1977, Hoyng, Knight and Spicer 1978). If we assume that 
this rate (2) is obtained after a linear increase over 10s, a typical time- 
scale for the impulsive phase, then our specified j b (t) is: 

j b (t) - 5.10 9 t statamp cm“ 2 (3). 
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The electron beam is decelerated by the electric field that drives the 
return current and ignoring electron time-of-f light effects, it has a "stopping 
length", s (Brown and Hayward 1981) due to ohmic losses, given by 


s 



(4) 


where E Q is the injected electron energy, here taken to be 100 keV and n is 
the resistivity of the plasma. Finally, we also assume a homogeneous and 
initially isothermal plasma with density n ■ 10 11 cm -3 and temperature 
T q - 5.10® K . 


The heating equations for the plasma electrons and ions, neglecting 
collisional heating by the beam, and also convection, thermal conduction and 
radiation losses, are 


3 

2 


dT e 

V dT 
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(Tj“T e ) 
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+ l Xihi.jp 
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( 6 ) 


where ic is Boltzmann's constant, t is the ion-electron energy exchange time, 

* X^CTg/T^) and (1-xp are the fractions of the ohmic power dissipation 
r^j 2 (t) absorbed by the electrons and ions respectively as a result of 
collision process i, and is the resistivity due to this process; e.g. in 
the case of classical Coulomb collisions, x s 1 and n is the usual Spitzer 
(1962) resistivity, proportional to T“3£. 

With the onset of turbulence, anomalous resistivity sets in due to wave- 
particle collisions and so the ohmic heating term £ n^j^ 2 rises, leading to 

i 

increased plasma heating. Note that this model, with a prescribed j (t), is 
fundamentally different from the Duijveman et al. (1981) problem where the 

electric field £(t) is specified. Specifically, their ohmic heating term, 
c 2 

— , is reduced by any turbulent increase in n, a fact they appear to have 
i 1 

overlooked. In other words, having a prescribed 5 -field leads to reduced, 
rather than increased, plasma heating when anomalous resistivity sets in. 

To solve equations (5) and (6) for T g (t), T^(t) and n(t) we use the 
hypothesis of marginal stability (Manheimer and Boris 1977) . States of 
marginal stability (i.e. states with zero growth-rate) in the present case can 
be represented by critical drift velocity (v ^ ) curves in the 




plane, where v 


- (^r ^ 


is the electron thermal speed. We apply 


v . at onset of turbulence, 
crxt * 


the marginal stability hypothesis by setting 

thereby constraining the system to evolve thereafter along the marginal 
stability curve. This gives us a relationship between T e and which, using 
the forms of x for ion-cyclotron (IC) and ion-acoustic (IA) waves derived by 
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Duijvetnan et al. (1981), then allows us to solve numerically the heating 
equations (5) and (6) . 


zero d-i] 

(a.i) 

\v 0 T ,/ 


We first consider the simplest case of an unmagnetised plasma, i.e, 
zero B-field. Figure 1 shows the evolution of the plasma in the 

plane. In the initial phase of classical resistivity, collisions 

'e *i' 

cause an increase in T e /T£, from its initial value of 1.0, while the drift 

velocity rises, until the IA marginal stability curve is reached. The system 

then evolves along the curve. Figure 2 shows the variation of the normalised 

resistivity with time, ion-acoustic turbulence switching on after about 6s. 

Finally, the increases in T e and Tj are seen in Figure 3. These profiles 

relate only to a plasma layer of thickness s . , 5.10 s cm (emission measure 

min 

n ^ s m £ n — S.lO^cm 3 ), the smallest depth of beam penetration into the 

atmosphere, because it is this region which is heated continuously during the 
simulation. The final temperatures obtained (T % 300 T 0 > T- 35 T ) will 
tend to be reduced by losses (neglected here) but we should note that wave 
generation will reduce thermal conductivity by the same factor as the electrical 
conductivity. 



Figure 1. 
stability. 


Evolution of the 


plasma in the 



plane, 


IA denotes the ion-acoustic marginal stability 


assuming marginal 
curve . 
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Figure 2. Variation of resistivity, normalised to the initial classical 
(i.e. Spitzer) resistivity. IA turbulence is generated at about 6s. 
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Figure 3. Rise in electron and ion temperatures, normalised to the initial 
temperature = 5.10 6 K, showing the increases at onset of IA turbulence. 



We now attempt to repeat the analysis for a magnetised plasma. For 
Tg/Ti < 8, the critical drift velocity for the generation of ion-cyclotron 
waves has the lowest threshold. These preferentially heat ions and so tend 
to reduce T e /T i . However they quickly saturate and do not contribute greatly 
to the total resistivity. The system then evolves away from the IC marginal 
stability curve, with T e /T £ once again increasing due to the dominant effect 
of Coulomb collisions, until the IA marginal stability curve is reached. In 
this case, onset of IA turbulence occurs at lower T g /T£ and we often find 
unphysical behaviour in the form of negative resistivity. Investigating 
further, by assuming a simple analytic form for the IA marginal stability 
curve and by using the form of ■ X^O^/T^ derived by Duijveman et al. 

(1981), we find that n is negative if Tg/Tj < 4.5. In other words, the 
application of marginal stability leads to negative values of anomalous 
resistivity, if the onset of ion— acoustic turbulence occurs at T /T. < 4.5. 
We discuss this failure of marginal stability further in Section 5. 1 


3 Wave Growth Analysis 

. * n tlie Previous discussion we did not incorporate any equations describ- 
ing the growth of the ion-acoustic waves, which occurs on very short timescales 
in comparison to , the rise time of the return current. Here we adopt 

the following equation for the evolution of W, the wave energy, 

dW 

dF - YW (7) 

where y, the linear growth rate, is taken to be 10" 2 a> p . ( Wp . is the ion plasma 

frequency). This growth rate corresponds to a drift velocity just in excess 
of the marginally stable drift velocity (Stringer 1964). Provided that 

w 

n kT <<: 1 * (i* e * weak turbulence), we can relate the effective collision 
P e 

frequency, v , and the wave level using: 


u 


W 


V - - w UJ 1 

eff pe n kT 
P e 

(see, e.g. Hasegawa (1974)) where w pe is the electron plasma frequency. 

effective (i.e. total, including classical) resistivity is then obtained 
using: 

n eff - ,. e ^ 


0) 


(see, e.g. Papadopoulos 1977). 


pe 


( 8 ) 

The 

(9) 


We can, therefore solve equations (5) and (6) numerically in conjunction 
th equations (7), (8) and (9). The wave level grows from the thermal level, 
obtained by substituting the classical collision frequency in the left hand 
side of equation (8). There are several possible saturation mechanisms limit- 
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ing the growth of ion-acoustic waves (Hasegawa 1974) and the relevant process 
for solar flare conditions is by no means certain. Hasegawa gives three 
possible processes, the saturation level being lowest for the nonlinear two- 
wave interaction investigated by Tsytovich (1971). The effective collision 
frequency at saturation for this process is 


where c 

s 



v eff * l0 ‘ 2 


0) . 
pi 



is the ion sound speed. 


( 10 ) 


Our approach then, is to let classical heating proceed as rises until 
the critical drift velocity for onset of ion-acoustic turbulence is reached. 

The wave energy then grows from the thermal level in the classical regime 
until saturation occurs. During the ion— acoustic heating phase, v d^ v e 

decreases due to the rapid increase in T e , until the system drops below the 
IA marginal stability curve (see Figure 4). The waves are then allowed to 
decay, with the same e-folding time (1/y) as the wavegrowth, until the wave 
level has dropped to the thermal level, where, once again, classical heating 
only takes place. 

Initial investigations show a rise in anomalous resistivity to ^ 10 3 times 
the initial classical value (see Figure 5), while T^ and T^ rise to ^ 50 T Q and 

^ 8 T q respectively (see Figure 6). These values are obtained in a very short 
time: the rise, saturation and decay of the IA waves takes place in ^10~^s . 
Again we note that the final temperature values are attained in the beam layer 
of thickness s m £ n » here found to be only 3.5.10 3 cm (emission measure 

— S.S.lO^cm -3 ), because of the very high resistivity attained during IA 
turbulence. In our simulations, it is found that saturation occurs at 

y q 

n k t ^ 10 ” thus satisfying the requirement for weak turbulence. 

P e 

The above analysis has, of necessity, been a very simple approach to the 
problem. For example, a wave energy equation with a proper averaging over 
the wavenumber, k, of the wave spectrum W(k) would be more realistic. In 
addition, the controversial issue of the relevant saturation mechanism for 
solar flare parameters needs to be looked at in more detail. 
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Figure 4. Evolution of the plasma in the pl ane using the wave 

\ v e A i/ 

growth analysis. The initial value of v^ is chosen such that the IA curve 
is reached after just a few numerical timesteps. 


NORMALISED RESISTIVITY AGAINST TIME 



Figure 5. The variation of normalised resistivity with time, showing the 
growth, saturation and decay of IA resistivity in 10 -5 s . 
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3: 




In this section we briefly compare the thermal and non-thermal 
bremsstrahlung emitted at the instant when the beam has shrunk to its 
minimum length s . n . From equation (4) it can be seen that when n rises at 

onset of turbulence, the beam length is reduced by the same factor (over the 
short time-scale involved, the current density j^Ct) is virtually constant). 

In other words, the electron beam lifetime is reduced and, correspondingly, 
the non-thermal bremsstrahlung. The enhanced ohmic dissipation of the return 
current, however, causes rapid local heating of the plasma and hence increases 
its thermal bremsstrahlung. 

Below, we calculate the ratio of thermal to non-thermal radiation in the 

case of beam return current losses only, which as we will show, in some 

circumstances, dominate collisional losses from the beam when n is anomalous. 

For comparison, we also evaluate the ratio of the thermal bremsstrahlung in 

the anomalous case to the beam bremsstrahlung evaluated from the usual thick 

target formula which incorporates collisional losses only. The thermal 

emissivity from the volume V = A.s . is 

min 

wbere tbe number density of electrons per unit energy range, is given by 
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the Maxwellian distribution function and (e,E) is the bremsstrahlung 

cross-section differential in photon energy. Here, as a simple first approx' 
imation, we assume Kramer's cross-section: 


^-(e.E) = 

8 

where Q = — ar ^ in the usual notation, 
^o 3 b 

emissivity (Brown 1971) is 


Qo m e c2 

eE 

The thick-target non-thermal 


( 12 ) 



f°° 

3(E)v(e,E) dE 

£ 


(13) 


where (E) = ^6(E-E o ) (s -1 per unit electron energy) is the electron 

flux and v(e,E) is the number of photons of energy e emitted by an electron 
with initial energy E, given by: 


v(e,E) 



de 


v(E.) 

^W) 


dEj 


(14) 


In the case of beam losses due to ohmic dissipation of the return current only, 
the energy loss equation is 

= -e£v(E*) (15) 

dt 

where the electric field, is given by 




en 


(16). 


The ratio of thermal to non-thermal emissivities (in the case of return 
current losses only) is then: 



= ( 8 ^ 
Virm kT ) 

x e e 7 



(17). 


Equation (13) also applies to the case of beam Coulomb collisional losses 
only. The energy loss equation in this case is 


dE* 

dt 


Kn p v(E*) 

E* 


(18) 


where K = 2ir e 4 A in the usual notation. The ratio of thermal bremsstrahlung 
in the anomalous case to the non-thermal emissivity in the case of Coulomb 
collisional losses only is then: 
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(19) 



8 V VK 

^“e<T e y/ ^ 0 (E 0 -e) 


Applying equations (17) and (19) to the results of Section 2 (Marginal 
Stability) and Section 3 (Wave Growth Analysis) we find the following: 


(a) Marginal stability analysis (Figures 1-3) gives, at e = 20keV: 



( 20 ). 


Here the thermal emissivity is over one order of magnitude greater than the 
non-thermal (return current) emissivity, and substantially exceeds that from 
a collisional thick target beam. Thus, when one considers a beam injection 
rate and reduces the beam area, A, until the return current goes IA unstable, 
the beam length and bremsstrahlung are greatly reduced as usually assumed (e.g. 
Hoyng, Knight and Spicer 1978) but enhanced thermal bremsstrahlung from the 
rapidly heated plasma exceeds the thick target bremsstrahlung which is produced 
by the same beam over a large area A. 


(b) Using the wave growth analysis (Figures 4-6), when unphysical 
behaviour is found in a marginal stab il i ty treatment, gives at e = 20keV: 



( 21 ) 


i.e. the thermal and return current non-thermal emissivities are comparable 
but the thermal emissivity is negligible in comparison with the collisional 
non-thermal emissivity. This is due to the small size of the emitting volume 
defined by s .. . In comparing these results with those of (a) , we should note 

that the evolution time in each analysis is very different - 20s in the 
marginal stability case compared to 10 -1+ s in the wave growth case. However, 
we anticipate that over longer evolution times the wave growth analysis will 
result in a greater volume of material heated to 'u 10 5 * * 8 K and a greater 
thermal emissivity (and emission measure) . 


5 The Failure of Marginal Stability 

We return now to the problem of why the application of marginal stability 
leads in some cases to negative values of resistivity. There appears to be 

confusion in the literature as to whether or not there is a lower limit of 

T g /T^ below which the ion-acoustic instability cannot be excited and therefore 
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below which the marginal stability curve used by us, by Duijveman et al. 

(1981) and by Holman (1985) for v crit (IA) is irrelevant. For example, Kaplan, 

Pikel'ner and Tsytovich (1974) state quite explicitly (p.53) that "... it is 
necessary that T g > 51^" (for the growth of ion-acoustic waves) . If this is 
correct, then it would be inappropriate to use our v cr £ t and X" IA curves 

below 5, and so it would not be surprising that unphysical results 

arise in this region. Similar statements are to be found in Melrose (1985) 
and in Stix (1962) (p.214) where a precise lower limit of Tg/T^ 3.5 is 

set. However all of these statements appear to be based on the assumption 
that v d « v e . In contrast. Fried and Gould (1961) and Kadomtsev (1965) 
state that the ion-acoustic instability can arise in an isothermal (i.e. 

T e /Ti 'v, l.o) plasma if the drift velocity is high enough: Vj > v g . 

By examining the roots of the linearised dispersion relation describing 
longitudinal plasma oscillations for drifting Maxwellians, we have ourselves 
confirmed the latter case: the so-called "ion— acoustic" mode will become 

unstable in an isothermal electron— ion plasma, described by drifting 
Maxwellians, if the relative drift velocity > 1.34v e . We hope to justify 
this in a subsequent paper. We are, of course, still left with the problem of 
explaining the failure of marginal stability for T /T. < 4.5. It is possible 
that the form of the X -function (Tange and Ichimaru^ 1974) for ion-acoustic 
waves is the source of this difficulty and we intend to investigate this. 

6 Discussions and Conclusions 

We have presented a simple analysis of the problem of electrostatic wave 
generation by a beam-driven return current. The marginal stability approach 
fails below an electron-to-ion temperature ratio T e /T£ ^ 4.5. We believe that 

this is not due to any lower limit of T e /T£ necessary for the generation of 
unstable ion-acoustic waves: the ion-acoustic instability will arise in an 

isothermal electron— ion plasma provided v, > v . 

d ~ e 

Both the marginal stability and wave growth analyses presented here allow 
us to calculate the uniform plasma heating only in the volume bounded by s 

. . min* 

the minimum stopping length of the beam. A 2-D treatment of the combined 

spatial and temporal dependence is required to extend the problem to plasma 
heating outside this region. Future work should also incorporate a less 
idealised injected electron spectrum (e.g. power law). However, we believe 
that the essential result presented here of enhanced ohmic return current 
dissipation leading to rapid plasma heating to hard X-ray temperatures in 
milliseconds (i.e. typical spike burst durations) or less, will remain unaltered. 
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